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A measurement of primordial non-gaussianity will be of paramount importance to distinguish 
between different models of inflation. Cosmic microwave background (CMB) anisotropy observations 
have set unprecedented bounds on the non-gaussianity parameter /nl but the interesting regime 
/nl 1 i® beyond their reach. Brightness-temperature fluctuations in the 21-cm line during the 
dark ages (2 ~ 30 — 100) are a promising successor to CMB studies, giving access to a much 
larger number of modes. They are, however, intrinsically non-linear, which results in secondary 
non-gaussianities orders of magnitude larger than the sought-after primordial signal. In this paper 
we carefully compute the primary and secondary bispectra of 21-cm fluctuations on small scales. 

We use the flat-sky formalism, which greatly simplifies the analysis, while still being very accurate 
on small angular scales. We show that the secondary bispectrum is highly degenerate with the 
primordial one, and argue that even percent-level uncertainties in the amplitude of the former lead 
to a bias of order A/nl ~ 10. To tackle this problem we carry out a detailed Fisher analysis, 
marginalizing over the amplitudes of a few smooth redshift-dependent coefficients characterizing 
the secondary bispectrum. We find that the signal-to-noise ratio for a single redshift slice is reduced 
by a factor of ~ 5 in comparison to a case without secondary non-gaussianities. Setting aside 
foreground contamination, we forecast that a cosmic-variance-limited experiment observing 21-cm 
fluctuations over 30 < 2 < 100 with a 0.1-MHz bandwidth and 0.1-arcminute angular resolution 
could achieve a sensitivity of order /nl^’ ~ 0.03, ~ 0.04 and ~ 0.03. 


I. INTRODUCTION 

Increasingly precise cosmic microwave background 
(CMB) [Ilia and large-scale structure laa measure¬ 
ments have zeroed in on a rather simple model of the 
cosmos, requiring only a handful of parameters. In par¬ 
ticular, initial fluctuations seem to be mostly scalar and 
highly gaussian HHZ!. They are well described by a sim¬ 
ple power-law spectrum, whose slope is consistent with 
a single scalar field driving inflation while slowly rolling 
down a very flat potential Eli- Proposed experiments 
like EUCLID [10] and PRISM [11] will measure the in¬ 
flationary parameters even more precisely, and possibly 
extract additional quantities, such as the running of the 
scalar tilt [H] . This should further constrain the form of 
the inflaton potential during the quasi-de-Sitter phase. 

While single-held inflation has the merit of simplic¬ 
ity, a plethora of alternative models remain consistent 
with current data [IMI]- The main characteristic that 
differentiates them from the simplest inflationary sce¬ 
nario is that they can generate significant primordial 
non-gaussianities (PNGs). The simplest form of PNG is 
a non-vanishing three-point function for the primordial 
curvature perturbation C, parametrized by a dimension¬ 
less amplitude /nl ~ (C^)/(C^)^- Single-field inflation 
leads to a small three-point function, corresponding to 
/nl 10“^ [181119]. Alternative models typically gener¬ 
ate /nl ~ 1, as a result of interactions with other fields 
EOlllI], higher-derivative terms in the Lagrangian El- 
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I24j . or other mechanisms m- Measuring /nl ^ 1 is 
therefore a natural target for future experiments to start 
significantly constraining the physics of inflation [25]. 

The best constraints on /nl to date are obtained from 
CMB studies [55], and are consistent with zero, though 
with a large uncertainty, ^5 — 40 depending on the 
shape considered. CMB measurements are now cosmic- 
variance limited in temperature down to the photon dif¬ 
fusion scale corresponding to multipole i ~ 2000. The 
anticipated improvement in polarization measurements 
is expected to only marginally tighten the constraints on 
/nl- Reaching the /nl ~ 1 frontier will therefore most 
likely require other data sets. 

Fluctuations in the brightness temperature of the 21- 
cm line of neutral hydrogen have the potential to open a 
new window on the high-redshift universe El EH]- This 
observable can in principle allow us to probe a fantas¬ 
tic number of modes, largely surpassing those available 
from CMB observations alone. First, fluctuations are un¬ 
damped down to the baryon Jeans scale (with wavenum¬ 
ber k ~ 300 Mpc“^) three orders of magnitude smaller 
than the photon diffusion scale {k ^ 0.2 Mpc“^). In 
addition, whereas CMB anisotropies probe a single sur¬ 
face, a line such as the 21-cm transition makes it possible 
to observe the early universe in tomography, and to co¬ 
add the information from each independent redshift slice. 
While the 21-cm line can in principle be observed all the 
way to cosmological reionization, at 2 ^ 10, the signal is 
cleaner at higher redshifts 2 > 30, the dark ages preced¬ 
ing the formation of the first luminous objects. We focus 
on this redshift range in this paper. 

The technical challenges to observe high-redshift 21-cm 
fluctuations are daunting, and will most likely require a 
telescope array on the far side of the Moon EH] , as well as 
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foreground-removal of the Galactic synchrotron radiation 
to an exquisite accuracy. We will not tackle these prob¬ 
lems in the present work, but focus on another source of 
contamination: the intrinsic non-gaussian nature of 21- 
cm fluctuations, even for perfectly gaussian initial condi¬ 
tions. 

While the baryon-photon fluctuations are highly linear 
at the epoch of last-scattering z ~ 1100, the perturba¬ 
tions in the cold dark matter (CDM) and baryon fluids 
have significantly grown by z ^ 50. Even if they remain 
small enough that no bound structure has formed yet, 
gravitational growth leads to a non-linear dependence of 
the density field on initial conditions. In addition, the 
21-cm brightness temperature depends non-linearly on 
the local baryon density, velocity gradient, and tempera¬ 
ture. Both effects give rise to a non-vanishing three-point 
function for the 21-cm brightness temperature, orders of 
magnitude larger than that resulting from PNG. Unless 
treated appropriately, this can jeopardize the usefulness 
of 21-cm fluctuations to measure PNG. 

Two previous studies have partially addressed this is¬ 
sue. Ref. m computed the bispectrum of 21-cm fluc¬ 
tuations resulting from non-linear gravitational growth, 
but treated it approximately as a confusion noise rather 
than a bias. Ref. [28] computed all contributions of the 
secondary bispectrum, but did not account for it in their 
final forecasts. In addition they only computed the bis¬ 
pectrum for specific triangle configurations. These two 
groups moreover get significantly different final results. 

In this paper we compute the primary and secondary 
bispectra using the flat-sky formalism. This accurately 
reproduces the full-sky calculation with a much lower 
computational cost, and greatly simplifies the analysis. 
We show that the shapes of the primary and secondary 
bispectra overlap significantly. Unsubtracted, the sec¬ 
ondary bispectrum would lead to a bias A/nl 10^. 
Even percent-level residuals after subtraction would lead 
to a non-zero non gaussianity of order A/nl 10. This 
warrants a Fisher analysis, fitting simultaneously for the 
amplitude of PNG and for nuisance parameters charac¬ 
terizing the residual secondary bispectrum after a best- 
estimate is subtracted. For a single redshift slice, we find 
that the uncertainty in /nl after marginalizing over the 
nuisance parameters is increased by a factor of ~ 3 — 6 
in comparison to an ideal case without secondaries. Fi¬ 
nally, we optimally combine redshift slices accounting for 
the smoothness of the secondary bispectrum as a func¬ 
tion of redshift. Our forecasts for a cosmic-variance- 
limited experiment targeting 30 < z < 100 with a band¬ 
width of 0.1 MHz and angular resolution of 0.1 arcminute 
are: (Tfiocai ~ 0.03, cTfequii ^ 0.04, and Ufortho ~ 0.03. 

^NL ' /nl ■'NL 

For the same angular resolution but a bandwidth of 1 
MHz our forecast is crfiocai ^ 0.12, cr.equii ~ 0.39, and 

Jnl /nl 

CTfortho ~ 0.29. 

^NL 

The paper is structured as follows, in Section H, we 
briefly review the basic physics of the 21-cm transition 
and the flat-sky formalism. In Section HI, we compute 
the different sources of non-gaussianities, both primordial 


and secondary. In Section IV, we forecast the potential 
signal-to-noise ratio reachable for cosmic-variance limited 
experiments. We conclude in Section V. 


II. BASIC ASSUMPTIONS AND NOTATION 


A. 21-cm brightness temperature 


The spin temperature of neutral hydrogen Tg is defined 
as usual through the ratio of the abundance of atoms in 
the triplet state and in the singlet state: 


= 3 
no 


( 1 ) 


where T* = 0.068 K = 5.9 //eV is the energy difference 
between the two hyperfine levels. The abundances no 
and ni can be obtained to high accuracy in the steady- 
state approximation by equating the rate of upward and 
downward transitions: 


no(C'oi-f i?oi) — ni(Cio-I-Rio), (2) 


where the Cij are the collisional transition rates, pro¬ 
portional to the gas density and dependent on the gas 
temperature Tgas, and the Rij are the rates of radiative 
transitions mediated by CMB blackbody photons (specif¬ 
ically, Rio includes spontaneous and stimulated emission, 
and Roi accounts for absorption). Since we are concerned 
with the dark ages preceding the formation of the first 
luminous objects, we do not account for transitions re¬ 
sulting from inelastic scattering of Lyman-a photons (the 
Wouthuysen-Field effect [30H32] b The steady-state ap¬ 
proximation is very accurate as Cij + Rij ^ H at all 
times |33|. In the limit ^ Tgas,Tcmb, valid at all 
times, the spin temperature is then given by 


Ts = Tcmb 


C 
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Cio + Aio^^ 


{Tgas — Tcmh), (3) 


where Aio is the Einstein-A coefficient of the hyperfine 
transition. 

The brightness temperature (or more accurately, the 
brightness temperature contrast with respect to the 
CMB) resulting from the resonant interaction of CMB 
photons with the hyperfine transition is given by 

T^T = {Ts-T,^b){l-e-^), (4) 


where the superscript “loc” emphasizes that this is the 
local brightness temperature, at the location of the ab¬ 
sorbing gas. 

The optical depth r is a function of the local neutral 
hydrogen density tiho = nH(l — iCe) (where Xe is the free 
electron fraction) and the gradient drVr of the peculiar 
velocity along the direction of propagation [Mj : 

— ^ \3 _ ^10 _ /r\ 

^ “ 32^ Ts H{z) + (1 + z)drVr ’ ^ ^ 
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where A 21 ~ 21 cm is the transition wavelength. The 
second term in the denominator accounts for the per¬ 
turbation of the optical depth due to the perturbed ex¬ 
pansion rate along the direction of propagatior0 The 
21-cm transition is optically thin (r ^ 1) in the regime 
of interest, so that t{Ts — Tcmb)- The brightness 

temperature observed today is then just 


rpohs 

^21 


'Y'loc 

^21 

l + z 


= T 


Ts — Tcmb 
1 + z 


( 6 ) 


where Tc x (Tcmb ~ Tgas) is the rate at which Thomson 
scattering of CMB photons by free electrons heats up the 
gas. Since Tp oc the fluctuations of the gas tem¬ 

perature are coupled to those of the free-electron fraction 
5x^ ■ In principle this equation should be solved simulta¬ 
neously with the evolution of 6x^, obtained by perturbing 
the recombination rate [37]. We find that neglecting 6x^ 
leads to errors of order ~ 10% for the linear evolution 
and we shall set Sx^ —>■ 0 for simplicity. With this simpli¬ 
fication, the equation for to second order is 


B. Fluctuations of the 21-cm brightness 
temperature 


(1 - (5b -b + i^^rc(5rg„ = 0. (9) 

gas 


We define Sy = —(1 + z){drVr)/H{z). Up to terms of 
third order in the fluctuations, the observed brightness 
temperature _(we shall drop the superscript “obs”) takes 
the form |37j^ 

T2 I = T2i( 1 -b (5„ -b si) + (7b (5b + 7 t '^Tgaa) (1 + ^v) 
+Tbb Sl -b TbT ^bSTg^^ + Ttt '^Tgas’ C^) 

where Sb = Sub/ub is the fractional fluctuation of the 
baryon density and St^^^ is the fractional fluctuation of 
the gas temperature, which affect T 21 through the col¬ 
lision rates. This equation neglects fluctuations of the 
ionization fraction Xe ~ 10“"^ at the redshifts of interest, 
as they lead to negligible fluctuations of T 21 which is pro¬ 
portional to (1 — Xe). We compute the coefficients in the 
above equation as described in Ref. m- They ought to 
be used for detailed prediction when actual data is avail¬ 
able. For this study, however, we shall make simplifying 
assumptions regarding the gas temperature fluctuations 
in order to keep calculations tractable. We now describe 
our approximations. 

The evolution of the gas temperature can be obtained 
from the first law of thermodynamics. Neglecting fluctu¬ 
ations of the CMB temperature and the effect of gravi¬ 
tational potentials, the full non-linear equation is m 


We shall consider scales larger than the baryonic Jeans 
scale: fc <C fcj ~ 300 Mpc“^. On these scales baryons 
behave just like CDM, so their evolution equation does 
not depend on Tgas- Given Sb, we can therefore solve 
for the gas-temperature fluctuations. We decompose the 
baryon-density fluctuation into a piece linear in the initial 
conditions S^^^ and a quadratic piece Sj^^^ resulting from 
non-linear gravitational collapse. We can then solve for 
the linear and quadratic parts of Sxg,,^ '■ 



Tcmb 

Tgas 



S^2) ^ 

rji 1 gas 

J- gas 





( 10 ) 

(5(^^(ll) 


Our final approximation is to assume that S^^^ is uni¬ 
formly proportional to the scale factor a, i.e. S^^^ (x, a') = 
[a'/a)Sl^\x,a), independently of the position x, and 


similarly that S^^^ cx a^. We then solve Eqs. (10) and 
0 starting at z = 1000 with vanishing initial condi¬ 
tions. The mean free-electron fraction Xg required for 
Fc and mean gas temperature Tgas are obtained from 
HyRec |3H||33|. This allows us to obtain three coeffi¬ 
cients C'i(z), C 2 (z) and C^iz) such that 



2. 1 + (5t,. 

3^'’ 1 + Sb 

Tc 

^cmb ^gas 

T 



( 8 ) 


^ This term is often mislabeled as a “redshift-space distortion”. 
It indeed has the same form in the optically thin limit, but has 
a qualitatively different origin ESj. On the one hand, redshift- 
space distortions arise from the Jacobian of the transform from 
real to redshift space to which the observer has access. On the 
other hand the local velocity gradient modifies the local expan¬ 
sion rate and hence the optical depth (or escape probability) of 
21-cm photons, regardless of the relative velocity between the 
observer and the emitter. This effect is analogous to the per¬ 
turbation of the Lyman-a escape probability which is one of the 
sources of perturbed recombination [331136) . 

^ Our definition of Sv differs from that of Ref. [33 by a minus sign. 


= C2(^:)[^I^^(a:,2:)]^ + C'2(z)(5f^(x,z),(13) 
which we show in Fig. 

The assumption that S^^^ oc a and Sj^'^ oc is not quite 
correct. Indeed this assumes that baryons behave exactly 
like CDM. In reality, they start with different “initial” 
conditions at z « 1000, after they decouple from the pho¬ 
ton fluid shortly after cosmological recombination: their 
overdensity is typically significantly smaller than that of 
the CDM on sub-horizon scales, and their velocity field, 
though comparable to that of the CDM in magnitude, has 
a very different scale dependence (hence leading to the 
relative-velocity effect [40]). Baryons therefore take some 
time to “catch up” to the CDM, and their growth rate at 
early times differs from Sb oc a, and is scale-dependent. 
Given that Thomson scattering maintains Tgas = Tcmb at 
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FIG. 1. Coefficients of the approximate decomposition of 
the gas-temperature fluctuations as a quadratic function of 
baryon density fluctuations: STg^^{x,z) ~ Ci{z)3l^\x, z) + 
C 2 {z)[Sl^\x, z)]'^ + C 2 iz)S^f^\x, z). At high redshift, Comp¬ 
ton heating is efficient and maintains Tgaa = Fcmb, with neg¬ 
ligible fluctuations, so Ci ^ C 2 ~ C 2 ~ 0 . At low redshift, 
the gas decouples thermally from the CMB and starts cooling 
down adiabatically, asymptoting towards Tgaa oc which 

implies Ci « C 2 —>■ 2/3 and C 2 —>■ —1/9. 


z > 200 , regardless of the exact value of 5^, this should 
not be a major issue, but should be properly accounted 
for in a detailed analysis. 

With these caveats in mind, we substitute our approx¬ 
imation a = into equation 

0 and obtain the following simpler expression for the 
21 -cm brightness temperature fluctuations to second or¬ 
der, with which we shall work for the rest of this paper: 

<5r2i « T2 i(<5W + > 5 ( 2 ) + [d(i)]2) 

+ a(z)dW(l + e)) +;3(^)[5(1)]2 

The effective coefficients a, (3, and 7 are straightfor¬ 
wardly obtained from the coefficients of Eq. Q and 
Ci,C 2 ,C 2 and are shown in Figure 


C. Neglected sources of fluctuations 

The above analysis is only valid on subhorizon scales, 
and does not account for several relativistic effects. First, 
the gas is not at rest with respect to comoving observers. 
We have already accounted for the resulting perturbation 
to the local Hubble expansion rate due to the velocity 
gradient. In addition, a local velocity leads to (i) a dif¬ 
ference between the proper time in the baryon rest frame 
and the comoving frame, (ii) a dipolar anisotropy of the 
CMB intensity in the baryon rest frame, and (in) an ad¬ 
ditional redshifting of the observed frequency. Gravita¬ 
tional potentials also affect the observed brightness tem¬ 
perature through: (i) a time dilation, {ii) a perturbation 
to the local expansion rate, {Hi) the Sachs-Wolfe and 


FIG. 2. Coefficients of the approximate decomposition of the 
21-cm brightness temperature given in Eq. (14 1 , in mK. In 
solid black we plot T 21 , in blue dashed a = dT2i/d5^\ in 
red dotted P = ^d^T2ild5l, and in purple dot-dashed 7 = 
dT2i/dSi^\ 


integrated Sachs-Wolfe effect^ and {iv) lensing by in¬ 
tervening structure, as is familiar from CMB studies. All 
these relativistic corrections are rigorously accounted for 
using the relativistic Boltzmann equation in Ref. [33]. 
They lead to fluctuations on scales comparable to the 
horizon at the redshift of absorption, i.e. k < 10 “^ 
Mpc“^ |33|. We will neglect them in this study, which is 
justihed as we shall see that most of the signal-to-noise for 
PNGs comes from small scales, with k ^ 10“^ Mpc“^. 

Redshift-space distortions are an additional source of 
non-linear fluctuations. The observer has only access to 
the total redshift Zobs = Aobs/A 2 i — 1 , and will compute 
the angular power spectrum on slices of hxed Zobs- The 
observed redshift is the sum of the cosmological redshift 
z and the redshift due to the relative peculiar velocity v\\ 
along the line of sight: Zobs = z -I- ti||/c. The observed 
brightness temperature at wavelength Aobs is therefore 

T 2 r(Aobs,n) = ^f^J^, (15) 

J- “T -S^obs 

where the true redshift z = Zobs — W|| (z, n)/c depends im¬ 
plicitly on the unknown local velocity. The angular power 
spectrum at fixed Zobs therefore has additional non-linear 
terms [531 133 SI] • We shall not account for those in this 
study but they should of course be modeled accurately 
when actual data is available. 

Finally, Ref. [37] showed that the non-linear depen¬ 
dence of the 21 -cm fluctuation on the local baryon den¬ 
sity and temperature leads to enhanced large-scale fluctu¬ 
ations due to the relative velocity effect [40]. The mag¬ 
nitude of the enhanced fluctuations is 5 T 21 ~ /3A((5g), 
where /3 is the coefficient of quadratic terms in the 


® See Ref. 1411 for a discussion of the ISW effect for 21-cm surveys. 
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brightness-temperature fluctuations, and A(^g) is the 
large-scale fluctuation of small-scale power due to the 
relative velocity effect. These enhanced fluctuations are 
most important for scales fc < 0.1 Mpc“^, and we will not 
account for them in this study, where we focus mostly on 
smaller scales. The relative-velocity effect also leads to 
a suppression of the average small-scale power, but this 
takes place at scales k > 100 Mpc“^, which we do not 
consider. 


D. Flat-sky formalism 


1. Fourier transform 


We consider a small patch on the sky, across which we 
can assume that the line of sight fi is a constant direction. 
We then define the Fourier transform of the brightness 
temperature as 


ST{k) = / drd‘‘x±e * '^dT{rn,x±). 


(16) 


Assuming matter domination and that the baryons have 
caught up to the dark matter so that 6 b oc a, at linear or¬ 
der the peculiar velocity term is 6 y{k) = (fc • n^Stik). 
The linear terms of Eq. (14) therefore contribute a 
Fourier transform 


5r‘“(fe) = [a + T 2 i{k-hy]Sb{k). 


(17) 


The power spectrum of 21-cm fluctuations is therefore 
anisotropic: to lowest order, and defining fc|| = k ■ n, 

PsAk) = (a + T2i kl/eyPsAk). (18) 

Similarly, the bispectrum ^ 2 , fcs) is anisotropic, 

and depends on the orientation of the wavenumbers with 
respect to the line of sight. It is defined as usual through 

{6T2i{ki)ST2iik2)ST2iik3)) = (27r)^(5D(fci + k2 + fca) 

xBsT{ki,k2,k3). (19) 


2. Harmonic transform 

Since we focus on small angular scales, we adopt a 
flat-sky formalism [43l|44]. We assume that the 21-cm 
temperature is observed with a finite window function 
W in frequency. The observed temperature is therefore 
the convolution of the underlying temperature with W, 
which we shall denote by W * ST. We define the flat-sky 
harmonic transform 


of sight. In terms of the Fourier modes of ST, this gives 

/ d^k 

^e^^'‘iiW(kyST(k)(27rySn(rk^-£), 

. 

where W(ky is the Fourier transform of the window 
function and we have defined 

SAt) ^ ( 22 ) 

(27r)2 r2 

The function peaks at the origin, with value (5 d( 0) = 
/sky/TT, where /sky is the fraction of sky subtended by 
the survey. If has a characteristic width A£ ~ (/sky)”^^^ 
and integrates to unity. Finally, a convolution of Jd with 
itself gives back. 

The covariance of ST{£) at equal r is given by 
{ST{e)ST*{l')) = J d^kA‘^A^SBirk3_ - £)SArk± - t) 

/ dk\ I 

-^\W\ykyPsT (A||,fe^).(23) 

For i we may approximate fcj_ « £/r in 

the inner integral. Carrying out the outer integral, we 
arrive at 

{STie)ST*{e')) « {27Ty~SA£' - £)Ce, (24) 

where [33] 

^\W\^ih)PsT{k\AM- (25) 

We show the flat-sky power spectrum Ci computed with 
different widths of the window function and for several 
redshift slices in Fig. 



FIG. 3. Flat-sky power spectrum Ce in the limit of infinitely 
narrow window function, for redshifts (top to bottom) z = 50 
(blue), a = 100 (black) and 2 = 30 (red). We also show the 
Ci at redshift 2 = 50 for a Gaussian window function of width 
0.1 MHz (blue dashed) and width of 1 MHz (blue dotted). 


STir,e)= ( *ST){rn,xA. (20) 

Ja 

where n is the line of sight, assumed constant over the 
small survey area A, and a;j_ is perpendicular to the line 


Similarly, the three-point function of 5T{£.) defines our 
flat-sky bispectrum: 


{ST(A)ST{i2)ST{t3)) = {2A'^5A£i+£2+A)Bi^t,t,. 

(26) 
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with [ 33 ] 


1. Local 


/ dki\\dk2\\ ~ ~ ~ , 

^(fcl||)W^(fc2||)W^(-fcl|| - k2ll) 

Bst {ki\\,ei/r;k2\\,£2/r) , (27) 

where we have dropped the dependence on in the k- 
space bispectrum since it is fixed by the triangle condition 
given fei = k^n + ii/r and k 2 = k 2 \\n + £ 2 /r. Note that 
we do not use the Limber approximation and perform the 
full integrals over A:||’s. 

We now describe and compute the different contribu¬ 
tion to 


III. BISPECTRUM OF 21-CM FLUCTUATIONS 


The simplest form of PNG is of the local type, where 
the primordial potential $ is a local non-linear function 
of a gaussian field 4 >'- 

$(a;) = (j){x) + (0^(3;) - ■ (29) 

This implies a non-vanishing bispectrum for $, given to 
lowest order by 

fc 2 , fca) = [P^{ki)P^{k 2 ) + 2 perm.] .(30) 

This form of the bispectrum peaks in the squeezed con¬ 
figuration (fci <C A :2 ^ fcs and permutations). 

Local-type PNG typically arises in multi-field inflation 
models, such as the curvaton model or modulated reheat¬ 
ing [49] . 


The bispectrum gets contributions from primordial 
non-gaussianities, which we would like to extract from 
the data, but also from secondary non-gaussianities, aris¬ 
ing from the non-linear relation between the observable 
and the initial conditions, even if the latter are perfectly 
gaussian. 

We only consider multipoles i > 100, which correspond 
to wavenumbers k > 0.01 Mpc“^, as most of the signal- 
to-noise ratio for bispectrum measurements is expected 
to come from small-scale modes. We therefore neglect 
the contributions of relativistic terms to the bispectrum, 
in particular the ISW-lensing bispectrum, which is the 
dominant secondary bispectrum for CMB anisotropies 

gSHIHl- 


A. Primordial non-gaussianities 


The contribution of PNG to the bispectrum of 21-cm 
fluctuations can be obtained to lowest order by only con¬ 
sidering the linear terms in Eq. (141, and assuming that 


they are linearly related to the primordial curvature fluc¬ 
tuations. The Fourier transform of the linear terms is 


given in Eq. (17). We define M{k,z) = db{k, z)/^{k), 
where $ = (S/^C is Bardeen’s gravitational potential. 
The bispectrum of brightness-temperature fluctuations 
gets a contribution 


3 

BP"“(fei, fc2, fca) = n(« + T2iL^})M{ki) 


X B$(fci, /C2, fca) 


(28) 


Equilateral 


PNGs of the equilateral type arise when there are non¬ 
standard kinetic terms in the inflation Lagrangian, which 
are included in the so-called P{X) models of inflation 
[23] . concrete examples of which are k-infiation Elisn] 
and Dirac-Born-Infield inflation [221 El]. In these models 
the effective sound speed Cg can be very different from the 
speed of light (c = 1 ), and the non-gaussianity parameter 
is related to this departure via = —(35/108)(cJ^ — 

1 ) m- 

This shape peaks when the three modes cross the hori¬ 
zon at the same time, and hence fci ~ /c 2 ~ ^ 3 - A good 
template for it is [55], 


sr“(^i-fc2,fc3) = 6/^f'A| - 


1 


{kik2k^ 


(kiklklY " 


5 perm. 


2 perm. 


l(31) 


where is the normalization of the power spectrum of 

$: P<i,(fc) = 


3. Orthogonal 

The “orthogonal” shape of PNG was defined 

in Ref. [53] to be orthogonal to the equilat¬ 
eral shape for the scalar product = 

J2ki,k2,k3 Its 

form is 


from primordial non-gaussianities, where = {ki-n)/ki. 

We now review the different shapes of the initial po¬ 
tential bispectrum i?$(fci, fc 2 , ^ 3 ) that we will consider 
in this paper (see e. g. Ref. [5^ for a larger variety of 
shapes). 




-3 


(fcifc 2 ) 


4-r 


-f 2 perm. 


(fci/c2A:3) " 


(fcifcpf)" 3 " 


5 perm. 


..(32) 
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The models of Galileon inflation |S3] and ghost infla¬ 
tion [55] predict very high values of In. general, 

in terms of the Lagrangian for the Goldstone boson tt 
during inflation, both equilateral and orthogonal shapes 
arise from cubic kinetic interactions, and the /nls are 
linearly related to the coefficients of the and 7 r(d7r)^ 
terms [55] . 

It is interesting to also mention the folded form of non- 
gaussianity jSS] , where the shape of the bispectrum peaks 
at flattened (folded) triangles (fci = ^2 = ^ 3/2 and per¬ 
mutations). Initial conditions different from the stan¬ 
dard Bunch-Davies vacuum would give rise to this kind 
of PNG (57). It can be expressed as a combination of the 
two above, as - B”***®) / 2 . 


4 ■ Directional dependence 


In some models where inflation is driven by a gauge 
vector field |58| or in solid inflation [59] there is an addi¬ 
tional form of PNG, that induces an extra dependence in 
the angle between the ki vectors. In this case the bispec¬ 
trum can be decomposed in Legendre polynomials [60] . 
where each component would be 


Bj^\ki,k2,k3) 


AJ) 

JNL 


(/ci )P$ (/C2 )'Pj (cos 6»12) 


-I- 2 perm. , 


(33) 


where Vj is the Legendre polynomial of order J, and 
012 is the angle between fci and ^ 2 , whose cosine can be 
expressed as cos 012 = (fc| — fci — A:|)/ 2 /cifc 2 . We consider 
J = 1, 2 and 3. 


B. Secondary non-gaussianities 

1. Non-linear gravitational collapse 

The growth of overdensities by gravitational collapse 
is a fundamentally non-linear process, leading to a non¬ 
vanishing 3-point function, even when starting from per¬ 
fectly gaussian initial conditions. The resulting bispec¬ 
trum can be computed from second-order perturbation 
theory (see e.g. Ref. [6T]). The correlation of two linear 
perturbations with a second-order density perturbation 
or normalized velocity divergence (0 = — V • v/H) takes 
the form 


(5W(fei)(5«(fc2)<5(2Hfc3))' = 2P(fci,fc2)PiP2, (34) 
(<5(i)(fci),5«(fe2)0('Hfc3))' = 2G(fei,fe2)PiP2, (35) 

where (...)' is the three-point function divided by 
(27r)^(5D(fei -I- ^2 + ks), and Pi = Ps{ki) is the power 
spectrum of the linear overdensity. The mode-coupling 
kernels F{ki,k 2 ) and G{ki,k 2 ) are both of the form 


EIlISlQ 

K{ki,k2) = C 1 -I-C 2 ki-k2 


^ ^ ) +C 3 (fcrfc 2 )^- (36) 


For a GDM-only universe, ( 01 , 02 , 03 ) = (|, I, |) for F 
and (|, I, |) for G. In reality, however, baryons start 
clustering after recombination, while the GDM overden¬ 
sities have already been growing since their scales en¬ 
tered the horizon. Their density and velocity fields at 
recombination are therefore very different and the sub¬ 
sequent growth factor of matter fluctuations is therefore 
not just D{a) oc a. It moreover has a scale dependence, as 
baryons, though they start with effectively zero overden¬ 
sity at recombination {Sb «C 6 c on sub-horizon scales), 
have a velocity comparable to that of the GDM, but 
with a different scale dependence. Different wavenumbers 
therefore grow at slightly different rates. The coefficients 
Ci in Eq. (36 1 are therefore in reality weakly dependent 


on redshift and, perhaps to a lesser extent, on scale [ 66 ] . 
We shall ignore these complications here and take their 
standard values. 

Assuming 6 b = 6 and using 6 v{k) = y?9{k) and 
(5i^^(fc) = pL^6^^\k) [note that this last relation only holds 
for the first-order perturbations], the bispectrum of 21 - 
cm fluctuations due to gravitational collapse is straight¬ 
forwardly obtained from Eq. (14): 


fc2, fcs) = 2(a + r2i/i?)(a + T2iiil) 

X (7F(fei, fc2) + T2i(/ri + M2)"G(fci, fc2)) P 1 P 2 

-|-2 perm. (37) 


2. Non-linear relation between brightness temperature and 
baryon density 

The relationship between the 21-cm brightness temper¬ 
ature and the underlying density and velocity field is fun¬ 
damentally non-linear, due to (i) the non-linear depen¬ 
dence of the optical depth on the local peculiar velocity 
gradient (r oc 1/(1 — 6 y)), (ii) the non-linear dependence 
of the spin temperature on the baryon density and tem¬ 
perature, and (in) the non-linear dependence of the gas 
temperature on the baryon density. Therefore even for 
a perfectly gaussian underlying density field, this non¬ 
linear mapping leads to a non-vanishing bispectrum. 

This contribution to the bispectrum can be obtained 
from the following three-point functions: 

{6b{ki)6b{k2M]{k3)y = 2 P 1 P 2 , (38) 

{6b{ki)6bik2)[6b6y]{k3)y = (pf -\- fj,l)PiP2, (39) 

{6biki)6bik2)[6l]{k3)y = 2^l^4P,P2, (40) 


^ These coupling kernels are derived in the sub-horizon limit. Since 
we are mostly interested in small scales we shall not concern 
ourselves with subtle issues regarding the squeezed limit of the 
gravitational bispectrum on horizon scales |63ti65| . 







where the superscript (1) is implicit in all the fluctua¬ 
tions. Using Eq. (14), the explicit expression for the bis¬ 
pectrum arising from the non-linearity of 5 T2i as a tracer 
is then 


fca) — (a + T2i/ii)(a -I- 
X (2/3 -|- a{^\ + /if) -f 2T2i/tf/if) P 1 P 2 

+2 perm. (41) 


The total secondary bispectrum is obtained by sum¬ 
ming Eqs. (37) and (41). Note that the bispectrum aris¬ 
ing from Eq. (37) requires the kernels F and G to be 
non-zero, whereas the bispectrum from (41) does not. 


C. Numerical evaluation and comparison 


Inserting the Eourier-space primordial bispectra 
Eq. (28) and secondary bispectra Eqs. (37) and (41) into 
Eq. (27), we obtain the harmonic-space bispectra in the 
flat-sky limit. 


We show the total secondary bispectrum in Eig. 
along with the bispectra resulting from local, equilateral 
and orthogonal PNGs. As found by previous authors 
[28] . we find that the secondary bispectrum is typically 
at least two orders of magnitude larger than the bispec¬ 
trum due to PNGs for /nl = 1- This order-of-magnitude 
difference can be understood quite simply: the ratio of 
secondary to primary bispectra is of order 


((5(5(5(5) S{z) 

^ ((5M/nl$) "" /nl4>' 




FIG. 4. Bispectra of 21-cm brightness-temperature fluctua¬ 
tions resulting from secondary non-gaussianities and different 
shapes of primordial non-gaussianity, with /nl = 1, at z = 50. 
The top panel shows the bispectra for equilateral triangles 
{JL = l\ = I 2 = £ 3 ). The bottom panel shows the bispectra for 
squeezed triangles {£ = ii = I 2 ^ £3 = 1/50). In dashed blue 
we plot local, in dotted orange equilateral and in dash-dotted 
green orthogonal non-gaussianity. In solid black we plot the 
secondary bispectrum. The bispectra are computed in the 
flat-sky approximation for an infinitesimally narrow redshift 
slice. 


We know that S{z = 0) ^ 1 at the non-linear scale 
k^hiz = 0) « 0.1 Mpc“^. Scaling back to z = 100 gives 
S{z = 100) ^ 10“^ at fc ~ 0.1, with an amplitude increas¬ 
ing logarithmically with wavenumber. The primordial 
gravitational potential is nearly scale-independent and of 
order $ ~ 3 x 10"^ We therefore obtain s^ec/^prim _ 
few X 100 for /nl = 1, consistent with our more detailed 
calculation. 

Note that our estimates for both the primary and the 
secondary bispectra neglected higher-order terms, for ex¬ 
ample terms of order ((5^(5^(5^) in the secondary bispec¬ 
trum. These terms are suppressed by an additional factor 
of order <5^ ^ 10“^ at 2 ; « 100, and are therefore compa¬ 
rable to the primary bispectrum only if /nl few times 
10“^. We will not consider them in this study, but they 
should be accounted for in a final data analysis aiming 
for a few percent uncertainty in /nl- 

In practice, we carry out the integrals up to some max¬ 
imum multipole /max corresponding to the resolution of 
the observations. 


IV. FISHER ANALYSIS 

A. Bias due to secondary non-gaussianities 

Assuming a single type of primordial non-gaussianity 
with bispectrum = /nl^^^^jj if secondary non- 

gaussianities were negligible the minimum-variance cubic 
estimator for /nl from a single redshift z would be [67l - 

IMl 

f MSI 

■'NL (|5prim /,prim^^ ’ 

where we defined 

( 44 ) 
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and the scalar product ( , )2 is constructed as 


(B\B^), = 47r/sky 




d‘^£id‘^i2d‘^i3 
(27r)4 


^(■^l + -^2 + f-s) 


jDi t)3 

/^tot/^tot/^tot ' 
^€2 ^^3 


(45) 


defined in Eq. (46 ^ The equilateral and orthogonal- 
type bispectra have more complex shapes in k and i- 
space than the local type, which is why their overlap 
with the secondary bispectrum decreases with increasing 
^max while it increases for the latter. 

In the next section we describe how to deal with these 
degeneracies. 


In this equation (7^°* is the total variance of ST{i), due to 
cosmic variance and all other sources of noise, including 
the instrument and foregrounds. We assume throughout 
that the noise can be approximately computed neglecting 
non-gaussian contributions to ST{£). In practice, all our 
results will be quoted in the cosmic-variance limit, i.e. for 
= Ce given in Eq. ([2|) 

We saw in the previous Section that the bispectrum 
resulting from secondary non-gaussianities is much larger 
than the one arising from PNGs, typically by two to three 
orders of magnitude for /nl = 1- Using the estimator, 
Eq. (43), would therefore lead to a bias 


(^prim, 5 sec)^ 

A/nL = = C, 


^^prim ^prim^ '^prim,sec'W ^prini^ ’ 


(46) 

where Cprim.sec S [—1,1] quantifies the shape overlap or 
degeneracy of the primordial and secondary bispectra [ge¬ 
ometrically, Cprini,sec U the cosine of the angle between the 
two bispectra, for the scalar product (45)]. The shapes 
of primordial and secondary non-gaussianity being dif¬ 
ferent, one may hope that their overlap is small [27^ 
However, assuming a cosmic-variance-limited experiment 
with an infinitely narrow window function and a resolu¬ 
tion of ^ 0.1' (corresponding to fmax = 10^), we find that 


GoCjSec 0.89, 


^equijSec 0.T9, and CQ,.tho,sec 0.83. 


The unsubtracted secondary bispectrum would therefore 
lead to large biases = 870, = 3900, and 

A/§l*'° = —3900. For maximum resolution of 1' (corre¬ 
sponding now to .^max = 10'^), the values of the degener¬ 
acy coefficients would be Cioc.sec = 0.80, Cequi.sec = 0.89, 
and Cortho.sec = —0.88, which in turn would make the bi¬ 
ases = 420, A/X‘ = 2400, and = -2100. 

Such a strong degeneracy may seem surprising at first, 
given the large number of triangles on which the scalar 
product depends. However, because the bispectra are es¬ 
sentially smooth featureless functions of £ for small angu¬ 
lar scales, they can have significant overlap in the sense 


B. Estimators for a single redshift slice 

One could in principle try and model the secondary 
bispectrum from first principles and subtract the result¬ 
ing bias A/nl from the estimated PNG amplitude. This 
strategy is the one adopted for the bispectrum of GMB 
anisotropies, where the main contaminant is the ISW- 
lensing bispectrum. Given the now well-measured cos¬ 
mological parameters, the latter can indeed be modeled 
to sufficient accuracy, i.e. with an error smaller than the 
statistical uncertainty in /nl [lamiTniin]- In the case 
of 21-cm fluctuations, however, even percent-level resid¬ 
uals in the modeled secondary bispectrum would lead 
to biases of order A/nl 10 , significantly larger than 
the statistical errors one may hope to achieve. Reaching 
sub-percent accuracy would require, first, a very careful 
treatment of subtle microphysical processes affecting the 
population of the hyperfine levels [72] • In addition, it will 
be limited by the accuracy of cosmological parameters. 

The amplitudes of all the secondary bispectra depend 
on the coefficients in Eq. ( |36| for the kernels F and 
G. Although we use the values derived for a matter- 
dominated, GDM-only universe in our analysis, we as¬ 
sume that the exact values could be computed exactly 
should one want to do so (Ml- On the other hand, the 


four parameters Ai = T 2 \,a, (3,^ in Eq. (14) depend on 


the detailed microphysics of the hyperfine transition. For 
now we assume that they can be modeled up to subper¬ 
cent accuracy and denote their best estimates by A^. 
Our model for the bispectrum is 

4 

AA = BTXl, + /nl + E /* > (47) 

i=l 

where is the best estimate for the secondary bis¬ 

pectrum obtained with the A°, fi = AAi are the un¬ 
known residuals of the four coefficients and = 

dAi. To make the notation more compact 
we denote^ /o = /nl and and recall 


^ Note that in order to have tenth-of-arcminute resolution at red- 
shift z = 100 one would need a baseline D > 350 km. In order to 
reach cosmic variance limit at 2 = 50 and for resolution of one 
arcminute, the parameters of the interferometer would have to 
be really optimistic, with complete coverage /cover = 1, a base¬ 
line of order the diameter of the moon D = 3500 km, and a time 
of observation of 2 years. 

^ Ref. |27| treat the secondary non-gaussianity as a source of noise 
instead of a bias, which is inappropriate. 


^ Consider for instance the 1-dimensional scalar product {F.G) = 
r^max with ^min ^ ^max- If F(i) oc and 

^min 

G{i) OC then their degeneracy coefficient is approximately 
c = \/l -h 2x/(l H-fc), where x = (^ —a)/ (2q;-|- 1). This is greater 
than 0.4 for x < 10. 

® Note that /nl and fi do not have the same dimensions, but this 
does not affect the analysis. 
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that we search for one type of PNG at a time. No¬ 
tice that we disregard higher order correction terms, 
proportional to and since we will be able 

to model those two parameters to a precision better 
than 0.3%, as shown in Tab. [T^ which would mean that 
the bias in the non-gaussianity amplitude is of order 
A/nl ~ 10^(3 X 10-3)2 < 10-2. 

We fit simultaneously for the amplitude of the PNG 
and for the residual coefficients of the secondary bispec¬ 
trum. We treat the latter as nuisance parameters over 
which we will marginalize. In geometrical terms, we con¬ 
struct an estimator for the PNG by projecting the ob¬ 
served bispectrum on the component of the primordial 
bispectrum orthogonal to all shapes of secondary non- 
gaussianity. 

The minimum-variance cubic estimators for the pa¬ 
rameters fi are given by m 

/, = (48) 

3 

where F~^ is the inverse of the Fisher matrix F whose 
components are 

F,, = {b\b>),. (49) 

The variances of these estimators are m 

<^1 = (50) 

and the signal-to-noise ratio (SNR) for /nl is 

/nl/\/ (.F-3)oO- 

We show in Fig.j^the forecasted SNR for the local-type 
PNG, for a single narrow redshift slice around z = 50, 
as a function of the maximum multipole moment £max 
(with £min = 100). We also show for reference the SNR 
one would obtain if one neglected the secondary non- 
gaussianities, i.e. when substituting {F~^)oo —>■ 1 /Fbo as 
in Ref. [^. We see that properly accounting for sec¬ 
ondary non-gaussianities and their correlation with the 
primordial bispectrum reduces the SNR by a factor of 
~ 6 . 

We also show the SNR integrated starting from = 
10 ® down to a minimum fminj as a function of the latter. 
It plateaus for t'min ~ lO^, so modes with smaller i do 
not contribute significantly to the signal-to-noise ratio, 
which justifies our neglect of several contributions to the 
bispectrum on large scales. 

In Fig. 1^ we show the forecasted SNRs for the 
other shapes of PNG we considered. Secondary non- 
gaussianities are less correlated with these shapes than 
the they are with the local type, so the reduction in SNR 
is not as dramatic (a factor of ^ 3). 

We summarize the forecasted SNR in Table Ufor a 
single narrow redshift-slice at z = 50, for irna,x — 10^ 
(corresponding of an angular resolution of roughly 1 ar- 
cmin) and £max = 10 ® ( 0.1 arcmin angular resolution), 
assuming a cosmic-variance-limited experiment (i.e. tak¬ 
ing (7)°* = Cl, and neglecting additional thermal noise). 



FIG. 5. Signal-to-noise ratio (SNR) for PNG of the local type 
with /nl = 1, for a single narrow redshift slice at z = 50 and 
assuming /s^y = 1. The blue dashed curve shows (Foo)'^'^^^ 
the SNR obtained if one neglected secondary non-gaussianity. 
The black solid and red dotted curves show [(F-^)oo]-^'^3, 
the SNR after marginalization over the unknown residual am¬ 
plitudes of the secondary bispectrum, as a function of 
(black solid) and as a function of ^min at fixed fmax = 10® 
(red dotted). 



FIG. 6. SNR for different shapes of PNGs (with /nl = 1), 
after marginalization over the residual amplitudes of the 
secondary bispectrum. The different lines correspond to 
equilateral-type PNG (solid black), orthogonal-type PNG 
(blue dashed), and the three direction-dependent shapes J = 
1,2 and 3 in dotted green, dash-dotted brown, and long- 
dashed red, respectively. 

In particular, we find that values of /^l 1-3 and ~ 0.23 

could be reached for .^max = 10 ^ and 10 ®, respectively. 
The bigger improvement for better resolution for the or¬ 
thogonal and equilateral shapes with respect to the local 
one is due to the fact that they become less degenerate 
with the secondary bispectra as more modes are added 
in the analysis, as argued in Section [IV A| 


It is interesting to discuss how well we could probe the 
four secondary coefficients, T 21 , a, /?, and 7 . In Tab.[ll|we 
show the relative errors reachable for each of them with 
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PNG type 

(arcmin) 

(t/jjl (0.1 arcmin) 

Local 

1.3 

0.23 

Equilateral 

14 

0.71 

Orthogonal 

11 

0.71 

J = 1 

83 

5.3 

J = 2 

4.5 

0.83 

J = 3 

40 

3.1 


TABLE I. Detection forecasts for different shapes of PNG for 
a cosmic-variance-limited experiment observing the full sky 
at a single narrow redshift slice at ^ = 50. The central col¬ 
umn gives the results for fmax = 10^ (equivalent to having an 
experiment with arcminute resolution), and the right column 
those for ^max = 10® (one tenth of arcminute). 


a 0.1 arcminute resolution, as well as the correlation with 
the rest of parameters. 



Tai 

a 

P 

7 

/nl 

Tai 

5.9 X 10“'‘ 





a 

-0.95 

2.6 X 10“® 




p 

-0.97 

0.91 

0.012 



7 

0.41 

-0.63 

-0.43 

7.2 X 10“'‘ 


/nl 

0.89 

0.85 

-0.92 

0.36 

0.23 


TABLE 11. Fractional error and correlation coefficients of /nl 
and secondary bispectrum amplitudes. The diagonal elements 
ar e the fr actional errors for each p aramete r, calculated as 
\J (^“ 1)00 for the case of /nl and sj (F“i)ii/Ai for the rest. 
The off-diagonal elements are the correlations between param¬ 
eters, defined as {F~^)ijl ^For these results 
we considered local non gaussianity, at redshift z = 50 and a 
resolution of 0.1 arcminutes. 


C. Choice of nuisance parameters 


In our analysis we have marginalized over the residuals 
of the four coefficients r2i,a,/3,7. Here we discuss how 
different choices would affect our results. 

On the optimistic side, if we were able to relate the 
four secondary coefficients to each other to high preci¬ 
sion we could choose to marginalize over a single overall 
amplitude for the secondary bispectrum. 

On the pessimistic side, we may choose to marginalize 
over all geometrically distinct contributions to the sec¬ 
ondary bispectrum. This would account for unknown 
redshift dependences in the Ci coefficients. Recalling 
that the kernels F and G are made of three geometri¬ 
cally distinct pieces, Eq. ( [W| ) gives 18 different geomet¬ 
ric shapes. Equation (411 adds three independent shapes. 


This amounts to a total of 21 distinct geometric shapes, 
the amplitudes of which we marginalize over. 

We show the resulting SNRs in Fig. where for refer¬ 
ence we also show the SNR in the absence of secondary 


non-gaussianities, and our main result, which considers 
4 nuisance parameters. As expected, our result lies be¬ 
tween the optimistic and pessimistic cases, which act as 
bounds for the SNR when considering additional sec¬ 
ondary bispectra. 

In particular, in the optimistic approach, the SNR is 
improved by a factor of ~ 5: we find detection limits 
^ 0.12, ~ 0.75, ^ 0.58, ~ 5.7, 

~ 0.88, ~ 13 at arcminute resolution and 

^ 0.0063, /^f0.032, /grtho ^ 0.030, ^ 

0.19, ~ 0.04, O-'IO at maximum resolution 

for a single redshift slice at z = 50. 

In a real experiment, a y^-like test should be carried 
out to find out whether additional secondary bispectra 
to the four proposed here need to be considered. 



FIG. 7. SNR for local-type PNG with /nl = 1 as a function 
of £max, when neglecting secondary non-gaussianities (top, 
solid black curve), marginalizing over an overall amplitude of 
the secondary bispectrum (blue dashed), marginalizing over 
4 coefficients as we do in the main text (red, dashed), and 
marginalizing over the amplitudes of the 21 geometrically 
distinct shapes of secondary bispectra (bottom, green dash- 
dotted). 


D. Tomography 


So far we have been studying the bispectrum on a sin¬ 
gle redshift slice, which would correspond to observing 
the 21-cm line with a single frequency channel. However, 
one of the great advantages of the 21-cm line is that it 
enables us to coadd information from different redshifts. 

Before thinking of how to add different redshift shells 
we will study whether they contain the same or different 
information. Let us construct a measure of the correla¬ 
tion between two slices at a radial distance Ar from each 
other. We define the correlation length ^r{^) as the radial 
separation beyond which the cross-correlation between 
two redshift-slices is less than 1/2 the power spectrum: 

Ce[Ar = m] = ic4Ar = 0], 


(51) 
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where the cross-power spectrum C£[Ar] is obtained from 

C 4 Ar] = ll^Psr {k^^,£/r) (52) 



FIG. 8. Correlation length as a function of £, defined as 
the separation in frequency beyond which two redshift slices 
are correlated by less than 1/2. This curve was calculated 
with an infinitely narrow bandwidth at 2 = 50, and for each 
i the correlation would increase to match the value of the 
bandwidth if it is bigger than the in the plot. 

From the correlation length in radial comoving sepa¬ 
ration we obtain the characteristic correlation length 
in frequency space through 

6 (53) 

az ar 

« 1 MHz (( —-—) , (54) 

\i + z) VeoMpcy’ ^ ’ 

where t'o = 1.4 GHz is the rest-frame frequency of the 
21-cm transition. 

We show the function in Figure]^ For I < 100 
(corresponding to A: < fceq ~ 0.01 Mpc^), P(fc||,£/r) 
peaks at fc|| ^ fceq, independently of and the cross¬ 
correlation C^[Ar] has a characteristic length scale « 
0.3 MHz, independent of i. For i > 100, the function 
P(fc||,£/r) has a characteristic turnaround scale at fc|| ^ 
i/r, which leads to a correlation length oc l/£. 

In order to compare with previous results in the liter¬ 
ature HZlllHl we will assume bandwidths Ai/ of 1 MHz 
and 0.1 MHz. As argued above (Fig. most of the sig¬ 
nal comes from large-£ modes {£ > 1000), for which the 
correlation length < 0.1 MHz, so in both cases we 
may assume that different redshift slices are completely 
uncorrelated. An observation of 21-cm fluctuations be¬ 
tween 14 MHz {z = 100) and 45 MHz {z = 30) with fre¬ 
quency resolution Av would therefore have Nances « 30x 
1 MRtj/Av independent redshift slices. 

The simplest analysis would consist in finding the best- 
fit /nl for each redshift slice and coadd the estima¬ 
tors with inverse-variance weighting. This procedure is 
not optimal, however, as the secondary bispectrum (and 


by extension, the residual after subtraction of the best- 
estimate is a smooth function of redshift. The 

redshift dependence of the residuals fi — AA^ can there¬ 
fore be modeled by a linear combination of a few basis 
functions and depends on a few coefficients instead of 
Asiices independent amplitudes: 

^bases 

Mz) = ( 55 ) 

1=0 

Several choices of basis functions could be made. We 
found that in the redshift range 30—100 the coefficients 
Ai{z) could be fit to ~ 1, 0.1, and 0.01 percent accuracy 
with third, fifth, or seventh-order polynomials in log(z), 
respectively. We assume that this will also hold for the 
residuals. We therefore adopt Pj{z) = [log(2/50)]'’ for 
j = 0 to Abases = 3 or 7 as our basis set. Our full model 
for the redshift-dependent bispectrum is therefore 

= Bll’Xiz) + / nl 6 ?™ 3 ( z ) 

4 -(Vbases 1 

+ (56) 

i=l 3=0 

where {z) = Pj (z) x (z). 

We now fit simultaneously for /nl and 4 x Abases nui¬ 
sance parameters fij. Because we assume the redshift 
slices are uncorrelated (specifically, the noise is uncorre¬ 
lated in different slices, but the signal is not), the total 
scalar product between two bispectra is simply obtained 
by summing the single-redshift scalar product over red¬ 
shift slices: 

(6^6'")=5](6^6'")., (57) 

Z 

where n = (ij) is a generalized index, and (6”,6'”)2 is 
the scalar product of two bispectra at redshift 2 defined 
in Eq. ( [4^ . The usual Fisher analysis leads to = 

(■F'“^)oo, where the (1 -f 4Abases) x (1 -f 4Abases) Fisher 
matrix Fnm is now defined from the total scalar product 
( [F^ . We find that using seventh-order instead of third- 
order polynomials degrades the SNR by no more than 
~ 20 %. The final results we quote are obtained using 
third-order polynomials. 

Our final results are shown in TablejTTlj where we quote 
the minimum /nl detectable for /sky = 1 and ^max = 10® 
for two different bandwidths (Az/ =1 and 0.1 MHz). For 
/sky < 1 all the results scale as oc 

In summary, with a bandwidth of 1 MHz we could 
cross the /nl = 0(1) threshold, enabling us to rule out a 
big class of models of inflation if no PNG is detected. In¬ 
creasing the frequency resolution to 0.1 MHz the numbers 
improve to /nl ~ few 10“^, which would be close to the 
ultimate limit of the consistency relation (/nl ~ ~ 1), 

and hence should be present even in the simplest model 
of inflation. 
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PNG type 

f^/NL (1 MHz) 

c^/nl (0-1 MHz) 

Local 

0.12 

0.03 

Equilateral 

0.39 

0.04 

Orthogonal 

0.29 

0.03 

J = 1 

1.1 

0.1 

J = 2 

0.33 

0.05 

J = 3 

0.85 

0.09 
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TABLE III. Minimum /nl detectable integrating all redshift 
slices between z = 30 and z = 100 for /sky = 1. In the central 
column we show the result for a bandwidth of Av = 1 MHz 
and in the right column for Av = 0.1 MHz. 


V. CONCLUSIONS 


Now that the information from the CMB on non- 
gaussianity has been almost fully mined, it is time to 
consider other potential data sets. Intensity fluctuations 
in the 21-cm line during the dark ages offer a window into 
yet unexplored times and scales, and a promising future 
probe of PNGs. 

The technical challenges that need to be overcome be¬ 
fore the required experiments see the light of day are 
daunting. Because of atmospheric attenuation it would 
require an observatory on the Moon. Even then, care 
should be taken with intense Galactic foreground emis¬ 
sion. Nevertheless, this is not an impossible task. 

An additional issue is that the 21-cm signal is intrinsi¬ 
cally highly non-gaussian, due to non-linear gravitational 
growth, and the non-linear mapping between brightness 
temperature and the underlying density field. In this pa¬ 
per we have, for the first time, addressed this issue with 
a rigorous Fisher analysis approach, assuming cosmic- 
variance limited experiments with a finite angular and 
frequency resolution. We have shown that for a single 
redshift slice the secondary bispectrum is significantly 
degenerate with the primordial one, which results in a 
noticeable decrease of the forecasted signal-to-noise ratio 
(SNR) for PNGs. This contrasts with the results of pre¬ 
vious work, where this degeneracy was either neglected 
when forecasting the SNR [28] , or where it was claimed to 
be weak m- We then co-added the information of inde¬ 
pendent redshift slices while enforcing a smooth variation 
of the secondary bispectrum amplitudes with redshift. 

For a full-sky experiment with Ajz = O.I MHz and 0.1- 
arcminute resolution, we forecast a sensitivity ~ 

0.03, which would enable us to check the famous inflation¬ 
ary consistency relation. We also forecast cr^equii « 0.04, 
« 0.04, « 0.1, /J=2 « 0.05, and /J=3 « 0.09. 

Measurements of 21-cm fluctuations therefore have the 
potential to significantly improve upon cosmic-variance- 
limited CMB bounds. 
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